INCOMPRESSIBLE MAGNETOHYDRODYNAMIC LIMIT OF THE 
VLASOV-MAXWELL-BOLTZMANN EQUATIONS 



XIANPENG HU AND DEHUA WANG 

Abstract. The hydro-dynamic limit of the Vlasov-Maxwell-Boltzmann equations is con- 
sidered for weak solutions. Using relative entropy estimate about an absolute Maxwellian, 
an incompressible Electron-Magnetohydrodynamics-Fourier limit for solutions of the 
Vlasov-Maxwell-Blotzmann equations over any periodic spatial domain in R 3 is studied. 
It is shown that any properly scaled sequence of renormalized solutions of the Vlasov- 
Maxwell-Boltzmann equations has fluctuations that (in the weak L topology) converge 
to an infinitesimal Maxwellian with fluid variables that satisfy the incompressibility and 
Boussinesq relations. It is also shown that the limits of the velocity, the electric field, 
and the magnetic field are governed by a weak solution of an incompressible electron- 
magnetohydrodynamics system for all time. 



1. Introduction 

The hydrodynamic models such as the Euler or Navier-Stokes equations were first es- 
tablished by applying Newton's second law of motion to infinitesimal volume elements 
of the fluid under consideration; while the kinetic equations are the mathematical models 
used to describe the dilute particle gases at an intermediate scale between microscopic and 
macroscopic level with applications in a variety of sciences such as plasma, astrophysics, 
aerospace engineering, nuclear engineering, particle-fluid interactions, semiconductor tech- 
nology, social sciences, and biology. If the particles interact only through a repulsive 
conservative interparticle force with finite range, then at low enough densities this range 
will be much smaller than the interparticle spacing. In that regime, the evolution of the 
density of particles F = F(x,£,t) is governed by the classical Vlasov-Maxwell-Boltzmann 
equaitons (VMB) pi QH EE] : 

a 77 

-zr+Z-V x F + e(E + £xB)-VeF = Q(F,F), x G M 3 , £ G M 3 , t > 0, (1.1a) 
at 

1 dE , , , 

- V x B = -/jtoj, divB = 0, on M 3 x(0,oo), (1.1b) 

OB o 

-|Vxfi = 0, div£ = — , on M 3 x(0,oo), (1.1c) 
dt Vo 

p = e Fd£, j = e F£d£, on M 3 x(0,oo), (l.ld) 
Jm.3 J R 3 
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where the nonnegative function F(t,x,£) is the density of particles with velocity £ at time 
t and position x under the effect of the Lorentz force 

E + ^xB, 

E is the electric field, and B is the magnetic field. The function j is called the current 
density, while the function p is the charge density. The constant e is the charge of the 
electron. The constant c is the speed of light. The coefficients fiQ and rjo are the mag- 
netic permeability and the electric permittivity of the plasma in the vacuum (see [18] ) , 
satisfying fiorjoc 2 = 1. The collison operator Q(F,F) is defined as 

Q(F,F)= [ def dLub(Z-Z*,cu)(F'Fl-FF*), 

where the nonnegative function 6(£, oj) given for £ G M 3 and oj G S 2 (the unit sphere in 
is called the collision kernel, and 



F*=F(t,x,&), F , = F{t,x,?), Fi = F(t,x,C), 

with 

£' = £-(£-£*,^V, 
£t = £* + (£ - £*,w)w, 

yielding one convenient parametrization of the set of solutions to the law of elastic colli- 
sions: 

/£' + £ = £ + £*, ( . 

\i£'i 2 + i£:i 2 = iei 2 + i£*i 2 . 1 j 

The interpretation of £, £*, £^ is the following: are the velocities of two colliding 
molecules immediately before collision while are the velocities immediately after 

the collision. We will consider the initial value problem of system (jl.ip with the initial 
condition: 

(F,E,B)\ t=0 = (F°(x,0,E°(x),B°(x)) for x G R 3 , £ G R 3 . (1.3) 

On the macroscopic level, the incompressible Electron-Magnetohydrodynamics-Fourier 
equations describe the evolution of the velocity field u = u(i, x) of an idealized fluid over 
a given spatial domain in R 3 under the magnetic field B = B(t, x) and the electronic field 
E = E(t,x), and take the form (cf. [1]) 

9 t u + u-Vu - pAu + Vp- aeE = (V x B) x B, (1.4a) 

d t B + VxE = 0, j = V x B = eu, (1.4b) 

d t 9 + u-9 = kA9, V x (h + 6) = 0, (1.4c) 

divu = 0, divB = 0, (1.4d) 

with 

« = - 1 TT / Kfexpf-ffW 
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3(2?r)i 
ssure, ' 
also be 

(u,B,6)\ t=0 = (u (x),B (x),8 (x)), x£R 6 , (1.5) 



where p, 9, E, h denote the pressure, temperature, electric field, and density respectively. 
The initial value problem will also be considered for system (|1.4|) with the initial data: 
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where 

uo, B G {v G L 2 (M 3 ) : divw = in £>'} and O G L 2 (R 3 ). 

We call (u,p,B,E,9) a loeafc solution to (|1.4j) - f)1.5[) if (u,p,B,E) is a Leray's solution 
of the incompressible electron- magnetohydrodynamic equation (|1.4a.j) - (|1.4b|) under the 
constraints (|1.4dp with initial data in fjl .5[) . while 9 is a weak solution in the sense of 
distributions to (|1.4c|) with the initial data in fjl .5[) . 

The motivation of this paper is to find a scaling and verify mathematically the transi- 
tion from the microscopic model (jl.ip to the macroscopic model (jl.4p as some parameter 
vanishes. One of the main objectives is to connect the DiPerna-Lions theory of global 
renormalized solutions of the Boltzmann equation with the Leray theory of global weak so- 
lutions of the incompressible fluid equations in a periodic spatial domain T = [0, l] 3 C M 3 . 
More precisely, we consider the hydrodynamic limit of the Vlasov-Maxwell-Boltzmann 
equations for weak solutions in this paper. Using relative entropy estimate about an ab- 
solute Maxwellian, an incompressible Electron-Magnetohydrodynamics-Fourier limit for 
solutions of the Vlasov-Maxwell-Blotzmann equations over periodic spatial domains in M 3 
is studied. It is shown that any properly scaled sequence of renormalized solutions of 
the Vlasov-Maxwell-Boltzmann equations has fluctuations that (in the weak L 1 topology) 
converge to an infinitesimal Maxwellian with fluid variables that satisfy the incompress- 
ibility and Boussinesq relations. It is also shown that the limits of the velocity, the electric 
field, and the magnetic field are governed by a weak solution of an incompressible electron- 
magnetohydro dynamics system for all time. 

We note that if the local conservation laws of momentum and energy are assumed, 
the similar result for the Navier-Stokes limit from the Boltzmann equation was verified 
in [I112]- Later, without making any nonlinear weak compactness hypothesis, Golse and 
Saint-Raymond in their breakthrough paper [16] established the incompressible Navier- 
Stokes limit of Boltzmann equations with bounded kernels, where they used the entropy 
dissipation rate to decompose the collision operator in a new way and used a new L 1 aver- 
aging theory to prove the compactness assumption. Recently, Levermore and Masmoudi 
|23j extended those results to a much wider class of collision kernels. As the collision 
is neglected, the Vlasov-Maxwell-Boltzmann equations become the Vlasov-Maxwell equa- 
tions. For Vlasov-Maxwell equations, the mathematical verification from the weak solu- 
tion of the Vlasov-Maxwell equations to the incompressible models has been conducted 
in [U [6j \TE[ 127] . When the solution is smooth, authors in [31 [22] considered the Hilbert 
expansion of solutions to the Vlasov-Maxwell-Boltzmann equations. For other related re- 
sults on the Boltzmann equation, see [131 El [23 EE] and the references therein. Also 
for the large-time behavior, stability, and regularity of solutions to the Vlasov-Maxwell- 
Boltzmann equations, see [91 \TU[ [201 [21] and the references therein. 

This paper is organized as follows. In Section 2, we state the formal scalings, the relative 
entropy, the technical assumptions, and the main result. Section 3 is devoted to a list of a 
priori estimates on the fluctuations of the density from the relative entropy. In Section 4, 
we consider the limit of the Maxwell equations. Section 5 will concentrate on the vanishing 
of conservation defects. And finally, in Section 6 we will give the detailed proof of our 
main Theorem 12.11 
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2. DlMENSIONLESS ANALYSIS, PRELIMINARY, AND MAIN RESULTS 

In this section, we first introduce the scaling of (jl.ip . then deduce the relative entropy, 
and finally state the main result. 

To begin with, we now focus on the nondimensional form of the Vlasov-Maxwell- 
Boltzmann equations. This form is motivated by the fact that the incompressible Electron- 
Magnetohydrodynamics-Fourier system (jl.4p can be formally derived from the Vlasov- 
Maxwell-Boltzmann equations through a scaling, when the density F is close to a spa- 
tially homogeneous Maxwellian M = M(£) that has the same total mass, momentum, and 
energy as the initial data. To this end, we introduce 



— &3~2 F ' E = T~E, B = — 



and 



-o, 



where the constants i*, x*, are the characteristic time, characteristic distance, and char- 
acteristic speed; see [5] for more physical interpretations of these constants. The non- 
relativistic effect requires 

= — , and e=(— J <C 1. 
t* \c J 

Substituting those new variables back to (jl.ip . and dropping hats, we obtain 
dF 1 

-K7 + i ■ V X F + e(E + e x B) ■ V^F = -Q(F, F), (2.1a) 
at e 

f)F 

e—-VxB = -j, (2.1b) 

^ + Vx£ = 0, (2.1c) 

div£ = 0, divE = (2. Id) 

where the coefficient e is usually refereed as the dimensionless mean free path or Knudsen 
number. 

Since the incompressible flow is the large-scale low-frequency fluid-like behavior of a 
plasma system ([51 H]), we need to further scale the time to the order of e" 1 . For this 
purpose, in the system (|2.ip . we further introduce the scaling as 

t = et, x = ex, £ = e£, 

F = \f, E=-E, B = -B, 

e° e e 

and 

b = e 2 b. 



Then substituting the above scaling back into (|2.ip . and dropping tildes, we obtain 
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dF 1 

£ ^7 + £ • + e£ ( eE + Z*B)- VtF = -Q(F, F), (2.2a) 
at e 

e^-V xB = - J -, (2.2b) 
dt e K ' 

^ + VxE = 0, (2.2c) 
div5 = 0, div£ = ^, (2.2d) 

The incompressible Electron-Magnetohydrodynamics-Fourier equations will be obtained 
when F is close to the absolute Maxwellian iW[jwith order e. Motivated by [2| 116} [23] . we 
set F = MG. Recasting the system (I2.2p for G yields 



dC 1 
e ~0^ + t ■ v xG + ee(eE + £ x B) ■ V^G - ee 2 E ■ £G = -Q(G, G), (2.4a) 

f) F i 

e^-V xB = - J -, (2.4b) 
dt e 

dB 

— + VxE = 0, (2.4c) 

divB = 0, div£ = £. (2.4d) 

where the collision operator is now given by 

Q(G,G) = [ [ {G'*G' -G*G)bfa-£,uj)duM,dZ*, 
Jr 3 Js 2 

where M* = M(£*). 

2.1. Relative Entropy. For any pair of measurable functions / and g defined a.e. on 
IR 3 x R 3 and satisfying / > and g > a.e., we use the following notation for the relative 
entropy 

which is a way to measure how far / is away from g. We are interested in the evolution of 

U £ {t) = eH{F £ \M) + - J (e\E £ \ 2 + \B £ \ 2 )dx, (2.6) 

where (F £ ,E £ , B £ )^ £>0 y are renormalized solutions (see definition in Section 2.4) of Vlasov- 
Maxwell-Boltzmann equations (|2,2p . This quantity contains the information from the 
standard (rescaled) L 2 norm of the electro-magnetic field and from the relative entropy 
between the renormalized solution F e (t,x,£) and the absolute Maxwellian M. 



H{f\g) = 11 
JtJm 



dfdxe[0,oo], (2.5) 



^The absolute Maxwellian is given as 

MiO^-^-r^pf-hl 2 ) (2-3) 

and corresponds to the spatially homogeneous fluid state with its density and temperature equal to 1, bulk 
velocity equal to and no effect from the electric field and the magnetic field. 
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The following lemma is devoted to the study of the evolution of the relative entropy, 
deduced from 

d 



dt 



H £ = e [ [ d t F £ {\nF £ - In M)d£dx 
JT JM 3 



(2.7) 



Lemma 2.1. Let (F £ , E £ , B e ) be a renormalized solution (refer to the definition in Section 
2.2.1 below) to (12.2ft . ThenH e (t) satisfies the differential inequality: 

Proof. In view of [21] , the inequality will follow from the lower semi-continuity of the weak 
convergence and an equality version when solutions are smooth. Thus, we will assume that 
those solutions are smooth. Observing that 

d t F £ In F £ = d t {F £ In F £ ) - 8 t F £ , 

we obtain 

e[ I d t F £ \nF £ dtdx = -^ f [ In (^^)(F £ i,F £ '-F £ *F £ )dtdx, 
JtJr 3 4e7 r y R3 \F £ *1< £ / 



d t F £ In Md£dx = -e 2 e I / F £ E £ ■ £d£dx. 
ITJM. 3 JT- 



and, by ([22 



Here, we used the following identity twice (see 

Q(/,/)C(£R = 7 / [ «* / duB(f'fi- f + -£}. 
4 JT il 3 JS 2 



IT 

Hence, 

' 77 1 ' 



/ / d t F £ (ln F £ - In M)d£dx = -±- [ [ In (F E '*F £ ' - F e .F e )d£dx 

JT JR 3 4e J T J R 3 \ F e ** £ J ^ ^ 



+ ee 2 / / F £ E £ -£,d£dx. 




it 



On the other hand, multiplying equation ()2.2b|) by E e , equation (|2.2cj) by B £ , integrating 
them in x over R 3 and then summing them together, we obtain, 

^ / (e\E £ \ 2 + \B £ \ 2 )dx = -- f E s ■ j £ dx = -e- I [ E £ ■ ^F £ d^dx. (2.10) 
dt J T e J T e J T J r3 

Substituting ()2.10|) back into (|2.9p . we obtain 

d t F £ {\nF £ - \nM)didx = -1 / / In (^f) {F^Fj - F £ *F £ )d£dx 
'TJR 3 4e JtJr 3 V F e* F e J 




e^d_ 
2 dt 



J^(e\E £ \ 2 + \B £ \ 2 )dx, 



which is exactly an equality version of (|2.8ft . □ 
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Notations. In order to avoid unnecessary constants in the sequel, we will assume that 
the nondimensionalization has the following normalizations: 

doj = 1, / Md£ = 1, 

S 2 Jr 3 

associated with the domain S 2 , and R 3 respectively; 

/ / G°Md£dx = 1, / / £G°MdZdx = 0, 

Jt Jr 3 Jt Jr 3 

[ f h\ 2 G°MdZdx= 3 -, 
Jt Jr 3 z z 
associated with the initial data; and 



I ff KZ*-Z,u;)du)M*(%*M(% = l, 

Jr 3 J R 3 J S 2 



associated with the Boltzmann kernel. 

Since Md^ is a positive unit measure on M 3 , we denote by < i] > the average over this 
measure of any integrable function r/ = ??(£), 



< 7] >= / r]Md^. 
Jr 3 

Since 

dM = b(C* - £,u))dwM*d£*Md£ 

is a non- negative unit measure on R 3 x R 3 x S 2 , we denote by C r > the average over 
this measure of any integrable function r = t(£, £*,cj), 



< r >= / rdM. 

Jr 3 

The collision measure dM is invariant under the transformations 

which are called collisional symmetries (cf. [2| I16j). 

Now, we can explain Lemma 12. II in terms of G £ as follows: 



d 



J {G £ \nG £ -G £ + l)dx+ E ^j t j (e\E £ \ 2 + \B £ \ 2 )d; 



£-77 J •<• in'" << • 1 • y-77 J H- |U-r ;c/./- 



(2.11) 



If G £ solves the VMB equations (|2.4p , then inequality (|2.1ip implies 

1 /•* 

H e (t) + - / K{G e (s))ds = U e (0), (2.12) 
e Jo 

where T~L £ (t) is the entropy functional 

H £ {t)=e [ <G £ lnG £ -G £ + l>dx+^ [ {e\E £ \ 2 + \B £ \ 2 )dx, (2.13) 

Jr 3 1 JR 3 
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and R(G) is the entropy dissipation rate functional 

n{G) =L\ ( in {Ge '* G£ ' ~ dx - (2 - i4) 

This choice of T~L £ as the entropy functional (|2.13[) is based on the fact that its integrand 
is a non-negative strictly convex function of G with a minimum value of zero at G = 1. 
Indeed for any G, 

H(G) > 0, and H(G) = if and only if G = 1. (2.15) 

Here H(G) is called the relative entropy with respect to the absolute equilibrium G = 1 
which provides a natural measure of the proximity of G to that equilibrium. 

We can expect that, the terms involving the entropy T~L £ measure the proximity of G e 
and G £ to the absolute equilibrium value of 1. On the other hand, the terms involving the 
dissipation rate 1Z, can be understood to measure the proximity of G £ to any Maxwellian 
through their characterization. 

2.2. Global Solutions. In order to mathematically justify the incompressible Electron- 
Magnetohydrodynamics- Fourier limit of the Vlasov-Maxwell-Boltzmann equations, we 
must make precise: 

• the notion of solutions for the Vlasov-Maxwell-Boltzmann equations; 

• the notion of solutions for the incompressible Electron-Magnetohydrodynamics- 
Fourier system (jl.4j) . 

Ideally, these solutions should be global while the bounds are physically natural. We 
therefore work in the setting of DiPerna-Lions renormalized solutions for the Vlasov- 
Maxwell-Boltzmann equations, and in the setting of Leray solutions for the incompress- 
ible Electron-Magnetohydrodynamics-Fourier system. These theories have the virtues of 
considering physically natural classes of initial data. 

2.2.1. Renormalized solutions to the Vlasov-Maxwell-Boltzmann equations. In the spirit 
of the DiPerna-Lions theory for the Boltzmann equation and the idea in Hu-Wang [21J, 
modified slightly for the periodic box, it is possible to show the weak stability of global weak 
solutions to a whole class of formally equivalent initial- value problems. More precisely, 
let G £ > be a sequence of DiPerna-Lions renormalized solutions to the scaled Vlasov- 
Maxwell-Boltzmann initial-value problem (j2.4[) with 

G £ (0,x,0 = G° £ (x,0>0, E £ (0,x) = E° £ (x), B e (0, x) = B° £ (x). 

A Renormalized Solution Relative to M of (12. 2p is a triplet (F £ , E £ , B £ ) such that 

F £ G C(R+;Ll oc (i? 3 ;L 1 (M 3 ))), E e , B £ G C W (R + ; L 2 (R 3 )), 

and satisfies 

r' Q(F £ , F £ ) G L} oc (R+ x R 3 x T) (2.16) 

for all r G C 1 (M + ) such that 



r(0) = 0, and z)T'(z) is bounded on R + , (2.17) 
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has finite relative entropy for all positive time: 

1 /"* 

U e {t) + - / K{G{s))ds < M e (0), (2.18) 
£ Jo 

and finally satisfies 



o JT 




DC 



JT 

oo 



3 



Ft 



e£ / / r 77 e£ £ ' CxMd^dxdt 



MJ >/v * (2.19) 



JT 




+ ee / / / r 77 ) £E e ■ CFeXd^dxdt 



M J 



1 -ri T L Y '{^) a(F " F ' )xMd(dxdt 



+ e* Jo 

= 

for each test function \ G Cq°((0, oo) x 

Throughout the rest of this paper, we assume that the initial data satisfies the 
normalizations and the entropy bound 



n e (0) < Ce 3 , (2.20) 

for some fixed C > 0. 



2.2.2. Weak formulation of the limiting system (|1.4p . Inspired by [29], for the limiting 
system (|1.4p with mean zero initial data, the Leray theory is set in the following Hilbert 
spaces of vector- and scalar- valued functions: 



H, = < 


Iw G L 2 (dx; 


M 3 ) 


K s = \ 


|x e L 2 (dx;] 


R): 


Vv = \ 


w G H„ : 




V H 


XGH S : 





divu) = 0, J wdx = ol , 



r 

Vwpdx < oo } , 



Let H = H„ © H„ © H s and V = V„ © V„ © V s . Leray's theory yields: given any 
(uo, Bo, 0o) G H, there exists a (u, B, 9) G C([0, oo); u> — H) fl Lf oc (0, oo; V) which equals 
initially (uq, Bq,9q) G H and satisfies the incompressible system (ll.4p in the sense that, 
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for all e HnC^T), 

/ • u(i)dx — / (/> • u(s)dx — / / V x : (u (8) u)dxdr 

Jt Jt Js Jt 

= -\i [ [ V x (f> : V x udxdr - ae [ [ Ecpdxdr + f [ V x (j) : (B ® B)dxdr; 

Js Jt Js Jt Js Jt 

f tp ■ B(t)dx - I ■ B(s)dx + I f E-(V x x ip)dxdr = 0; 
JT JT J s JT 

X0(t)dx - [ x9(s)dx - [ [ V x x • {u6)dxdr 
T JT Js JT 

/ V^x • V x 9dxdr, 
s Jt 

for every < s < t. Moreover, (u, B, 9) satisfies the dissipation inequalities 

j i(| u (t)| 2 + a|B(t)| 2 )+ f* n\V x u\ 2 dxds < I i(|u | 2 + Q|S | 2 )dx, (2.21a) 
[ -\6(t)\ 2 dx+ [ f n\V x e\ 2 dxds< [ -\6 \ 2 dx, (2.21b) 

Jt 2 Jo Jt Jt 2 

for every t > 0. 

A global existence theory, similar to Leray's theory of incompressible Navier-Stokes 
equations, can be established via Garlerkin's method, the dissipation inequalities (12.211) 
and Ohm's law which expresses the electric field E in terms of the magnetic field and the 
velocity as, see [51 [18] 

j = a (E + u x B), 

where a > is the electrical conductivity. To obtain the dissipation inequality (|2.21al) . 
we first multiply ()1.4aj) by u to obtain, using (|1.4bj) . 

1 d f 

2^ll u lli 2( r)+Mll Vu lli2 (r) -a J^E ■ (V x B)dx = 0. (2.22) 
Here, we used the identity 

B x (V x B) ■ u = -(B x (V x B)) - j = -(B x (V x B)) ■ (V x B) = 0, 

e e 

according to (|1.4b|) . Then, we multiply (|1.4bj) by aB to obtain 

%J t W B Wl\T) + a J r E-(VxB) = 0. (2.23) 

Adding (|2.22[) and ([2.23p . and then integrating it over (0,T) yield the energy inequality 
(|23Tal) . 

In summary, we have the following existence theory for the incompressible system (fTT 



Proposition 2.1. For each u ,B G {/ G L 2 (M 3 ) : div/ = in P'} and O G L 2 { 
there exists at least one weak solution (u,B,9) of (11 .4[> - (jl .5f) that satisfies the energy 
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inequality 



~ jf (ju(t, x)\ 2 + a\B(t, x)\ 2 + ||0(t, x)\ 2 ^j dx + jf_ (^/i|Vu| 2 + ^|V0| 2 ^ dxds 
i| (|uo| 2 + a|B o | 2 + ^|0o| 2 )^ 



< 
~ 2 



for all i > 0. 



2.3. Assumptions. In this subsection, we state our technical assumptions. To begin 
with, we define 



A(£)= / b{£,u)dw. 
Js 2 

Our assumptions regarding the collision kernel b are stated as follows: 

• (HO) b € L 1 {Br x S 2 ) for all R € (0, oo), where B R = {z G K 3 : \z\ < R}, and 

|6(z,w) depends only on |z| and |(z,o;)|, 

(^(l + l^l 2 )- 1 (j z+BR A(t)d£) -+0, as |*| -"►«), for all i?e(0,oo). 

• (HI) ^ < 6(z,w) < 6oo, zel 3 , oj £ S 2 , for some 6^ > 0; 

The assumption (HO) is assumed to make possible the global existence of renormalized 
solutions to the Vlasov-Maxwell-Boltzmann equations, see [8j[2l]. The class of collision 
kernels satisfying (HO), (HI) is not empty since it contains at least all collision kernels of 
the form b(z,oj) = b(\ cos(z,w)|) satisfying (HO). 

Next, we impose one more technical assumption on the sequence of fluctuations {ge}{ £ >0} 
(see (j2T25|) below). 

2 

• (H2) The family (1 + |£| 2 )-^r is relatively compact in w — L 1 (dtMd^dx), where 
N £ = l + lg £ . 

This assumption is the same as (A2) of Lions-Masmoudi |24j and similar to (H2) of [2], 
with the only difference being that we had to add the time variable, since we are dealing 
with the nonstationary case, when compared with the stationary case in [2]. 

2.4. Main Result. We consider a sequence of solutions G £ to the scaled Vlasov-Maxwell- 
Boltzmann equations 

ed t G £ + C ■ V X G £ + ee(eE £ + £ x B £ ) ■ V ? G £ - ee 2 E £ • £G £ = -Q(G £ , G £ ), (2.24) 

e 

in the form 

G £ = 1 + eg £ . (2.25) 

We expect that as e tends to zero, the leading behavior of the fluctuations g £ is for- 
mally consistent with the incompressible Electron-Magnetohydrodynamics-Fourier equa- 
tions. Indeed, formally, substituting (|2.25p into (|2,24p . we obtain 

edtge + £ • V x g e + ee(eE + (xB)- V^ e - eeE e ■ £ - ee 2 E ■ £g £ + -Lg £ = Q(g £ ,g £ ), (2.26) 

£ 

where L, the linearized collision operator, is given by 



Lg = -2Q(l,g)= f [ (g+g.-g'-g'^ 



)bduM*d^ 
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Repeated applications of the dA/f-symmetries yield the identity 
<vLg> = (v(g + g* - g' - g'*} 

= \ {{v + v*-v' - + g*-g' - <?*))} , 

for every v = i>(£) and g = g(£) for which the integral makes sense. This shows that L is 
formally self-adjoint and has a non-negative Hermitian form. Furthermore, using the dM- 
characterization, it can be shown that for any g = in the domain of L, the following 
statements are equivalent: 

Lg = 0; (2.27a) 

g = a + f3 ■ i + ]p\i\ 2 , for some (a, 0, 7) G R x R 3 x R. (2.27b) 

This characterizes N(L), the null space of L, as the set obtained by linearizing about 
(a,/3,7) = (0,0,0). From (|2.26p . we deduce formally that the limit of Lg e is zero and it 
can be expected that the limit of g £ will belong to N(L). Indeed, it was proved by Grad 
(see I16j) that for any collision kernel b satisfying (HI), L is a bounded nonnegative 
self-adjoint Fredholm operator on L 2 (Md^) with null space 

iferL = span{Ui,6,£ 3 ,l£| 2 }- 

Notice that since each entry of the tensor £ ® £ — ||£| 2 -f and of the vector ^£(|£| 2 — 5) is 
orthogonal to KerL, there exist a unique tensor $ and a unique vector such that 

L$ = £®£- i|£| 2 /, $ e {KerL) L C L 2 (Md£); ( 2 - 2 8) 

= -£(|£| 2 - 5), * G (KerL) 1 - C L 2 (Md£). (2.29) 
Now, our main result can be stated as follows. 

Theorem 2.1. Under the hypotheses (H0)-(H2), let G £ (t,x,£), with the form (|2.25|) . 6e 

a sequence of non-negative renormalized solutions to the scaled Vlasov-Maxwell-Boltzmann 
equations fj2.4j) satisfying the initial condition ()2.20p . Then, 

• T/ie sequence g e converges in the sense of distributions and almost everywhere to 
a function g as e tends to zero, and g is an infinitesimal Maxwellian, 

= ^ + u-£ + 0(~|£| 2 -~Y (2.30) 

where the velocity u satisfies the incompressibility relation, while the density and 
temperature functions, h and 9, satisfy the Boussinesq relation: 

divu = 0, V x (h + 6) = 0. (2.31) 

• As e — > 0, E £ and B £ converge to E and B in the sense of distributions and 
Lj°(L 2 (T)) respectively. 

• Moreover, the functions h, u, 9, B, and E are weak solutions of (jl.4p with 

1 2 
[i = — < $ : L$ >, k = — <^-1A>>. (2.32) 
10 15 
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3. Implications of the Entropy Inequality 

In this section, we first recall some results in [2j [16] which were established in the 
greatest possible generality, and relied only on the a priori estimates and in particular 
have nothing to do with the equations. To this end, from now on, we assume that the 
initial data G £ satisfies the entropy bound: 

e J <G°lnG°-G° + l>dx + y J {e\E° £ \ 2 + \B°\ 2 )dx < Ce i (3.1) 

with C > 0. From the relative entropy, we can obtain the uniform bound ||-Be||L°°(.L 2 (<£E))! 
and hence we can assume 

B £ ^ B weakly* in Lf (L 2 (dx)) , (3.2) 

with divi? = in T>' . Furthermore, from the relative entropy, H-Eell^^^)) is uni- 
formly bounded, and hence, we can assume that 

e^E £ ^n, weakly* in Lf(L 2 (dx)) (3.3) 
for some function £2 G L^°(L 2 (dx)). Then the results in [2J[TB], combining with (|2,11[) and 



(|3.1j) imply the following convergence. 

Theorem 3.1. Under assumptions (H0)-(H2), let F £ be a family of renormalized solu- 
tions to (|2.2p with initial data (F £ ,E £ ,B £ ) satisfying (|2.2(jp . and define the associated 
family of fluctuations by 

F £ -M 

Then 

• g £ is relatively compact in w — Lj oc (dtdx; L 1 ((l + |£| 2 )M<i£)) ; and for almost every 
t G [0, oo), g satisfies 



J \ (g 2 (t)) dx < liminf J /±h(eg E (t))\ dx < C; 



(3.4) 



moreover, for almost every (t,x), g(t,x,-) £ N(L), which means that g is of the 
form 

g(t, x, •) = h(t, x) + u(t, x) ■ £ + 6(t, x) (\\i\ 2 - I) , (3.5) 



2"' 2 



where (h,u,9) £ L°°(dt; L 2 (dx; R x R 3 : 
the rescaled collision integrands 

qe = ^(G' £ G' e *-G e G e *) (3.6) 

satisfy that ^{G £ )q £ is relatively compact in w — Lj oc (dtdx; L 1 ((l + |£| 2 )dM)); fur- 
thermore, any of the limit points q of r y(G £ )q £ as e — > satisfies the dM-symmetry 
relations 

= + <t>* - 4>' - #)<f) , (3-7) 

and, q 6 L 2 (0, T; L 2 (dMdx)) . 
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• for any subsequence e n — > such that 

g £n ~>g, and ~f(G £n )q £n q 

inw- L} oc {dtalx;L l {(l + |£| 2 )Md£)) and in w - Lj^dtdx; L 1 ((1 + \£\ 2 )MdM)) 
respectively. 

• denoting N e = | + \G e , then j^- is bounded in Lf 2 (L 2 (Md^dx)) and relatively 
compact in w - Lj oc (dtdx; + |£| 2 )dM)). 

Weak compactness statements regarding g £ and q £ result in the following bound for 
their limits. 

Lemma 3.1. Under the same conditions as Theorem \3.1[ for almost every t 6 [0, oo) the 
function g and q satisfy 

-t 



f I (/(f)) dx + 1 _( m 2 + \B\ 2 )dx + I f \ (q 2 » dxds 



<Uminfy /Lh(£g°)^ dx <C. 

Proof. Taking the lim inf on the both sides of the entropy inequality (|2.1ip , we obtain 

liminf / ( -^hUgJt)) ) dx + - liminf / (e\E £ \ 2 + \B e \ 2 )dx 
Jr\ £ I 2 Jr 



(3.8) 



+ lim inf 

£->0 



1 // 1 



o JT 



4 \\ e 4 V G 



2 



T *eG E 



G* F G F )) dxds 



(3.9) 



< liminf ^h(sG° £ (t))J dx + \ J^(e\E°\ 2 + \B° £ \ 2 )dx^j 



< C. 



Due to the lower semi-continuity of the weak convergence, we deduce that 

/ (\tt\ 2 + \B\ 2 )dx < liminf / (e\E £ \ 2 + \B £ \ 2 )dx, 
Jr Jr 

while, from the second assertion of Proposition 3.1 in [2], we have 
/ \{9\t))dx+ f f\(q 2 ))dxds 

liminf / / —= h(eq F (t)) \ dx + liminf / / 
^0 J T \£ 2 I e-+o Jo J 7 



(3.10) 



! ' ^ G* F G F \\ dxds. 



Substituting (j3TT0j) and (jBTTTT) back into (pT9j) . we finish the proof of (|3T8 



(3.11) 



□ 



To better understand the behavior of the fluctuation {g £ }{ £> o}, as in [IB] we introduce 
a class of bump functions 



T = | 7 :K+^ [0,l}\jeC\ 7 
We decompose g £ as 



3 5 
4' 4 



{1}, supp7 C 



1 3 
2' 2 



ff £ = 5* + ^ 



(3.12) 
(3.13) 
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with 

9 b e = - £ ( G e- lh(G e ), g c £ = ^(Ge- 1)(1 - 7(G £ )), 

where 7 G T. The following entropy controls (Proposition 2.1 and Proposition 2.7 in |16j) 
will be very useful: 

Lemma 3.2 (Entropy controls). Assume that the bump function 7 E T as in (|3.12p . 

TTie relative fluctuation g £ of the density satisfies the following estimates: 

g b £ = 0{l) in Lf{L 2 (Md£dx)); 

• (1 — j(G e )) < Ae 2 \g^\, which implies that ~(1 — 7 (G e )) < 2\g £ \z , and 

g £ = 0(l) in L™{L l {Mdidx))- 

• (1 - 7(G £ ))G £ < he 2 \g% and (1 - 7 (G e )) < 4e 2 | 5e c |. 



4. Implications of the Maxwell Equations 

For the asymptotic behavior of the solutions under the hypothesis H e (0) < Ce 3 , one 
of the difficulties when we deal with the magnetic field and the electric field comes from 
the fact that the relative entropy does not provide useful information on the electric field 
E £ due to the e in the front of the electric field in the definition of the relative entropy 
H. £ . Fortunately, the uniform estimate from the relative entropy is enough to ensure that 
Q = 0. Indeed, from p , 

e— - -> 0, in V'(R + x M 3 ). 
dt \ + j 

Next, since g £ converges to g in w — Lj oc (dtdx; L 1 ((1 + \£\ 2 )Md!;dx)), by the Cauchy- 
Schwarz inequality, we can deduce that g £ converges to g in w — Lj oc (dtdx; L (\£\Md£dx)). 
Due to the fact < £ >= 0, 2f =< ^g £ >, j| converges to j in w - Lj oc (dtdx). Then we 
take the limit as e — > in the equation (|2.2bj) to get 

VxB = j (4.1) 

in the sense of distributions. Furthermore, 



— II < £,9e >IIl°°(L 2 (T)) 



< l|ffs[U°? (L 2 (Md£dx)) (l£| 2 ) 2 

< OO. 

This implies that & converges weakly* to j in L£°(L 2 (T)). 
On the other hand, for the electric field E £ , we have 

Lemma 4.1. T/ie family {-E e }{ e >o} formally satisfies 

E £ = d t (eE £ x fl 6 ) -(VxB £ )x5 e + ^xB E + ediv(£ e ® £ e ) 

1 f £ (4.2) 

-e-V|£ £ | 2 -e£ E / g £ MdC, 

4 7lB3 
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in the sense of distributions. Hence, {E £ } £> q is uniformly bounded in (Wq ,oo ((0, T) xT))'. 

Proof. Indeed, multiplying (|2.2bp by B £ , multiplying (|2.2c[) by eE e , and adding them 
together to yield 

d t {eE £ x B £ ) - (V x E> £ ) x B £ + sE £ x (V x E £ ) =~x B £ . (4.3) 

Note that 

EdivE + (V x E) x E = d\w(E (£> E) — ^V\E\ 2 . (4.4) 

The identity (14. 3p can be rewritten as, using (I2.2d|) 
E £ p £ = eE £ divE £ 

U 1 o (4-5) 

= 9 t (e£ e x B £ ) - (V x B e ) x B £ + — x B £ + ediv(£ £ £ £ ) - e-V\E £ \ 2 . 



Because 

Pe 

one obtains, according to ([4.5 



/ (1 + eg £ )Mdi = [ Mdi + e [ g £ Md£ = 1 + e [ g £ MdC, 



E £ = d t {sE £ x - (V x B £ ) x B £ + j x B £ + sdw(E £ ® E £ 
- e\v\E £ \ 2 - eE £ [ g £ Mdi. 

2 ./W3 



(4.6) 



Next, due to the uniform bounds 

HV^^elU 00 (0,T;L 2 (R 3 )) < C, ||5e||L°°(0,T;L 2 (R 3 )) < C> 

we have 

in (W 1,oo ((0,r) x T))' as e -> 0, and -(V x B £ ) x B £ + J f x B £ is uniformly bounded in 
(W 1,oo ((0, T) x T))' by using the identity (g2D for 5. 
Also, we can control the term eE £ j R3 g £ Md£ as follows 



eE £ / g £ Mdi 



L 1 ((0,T)xT) 

l 

/" Jl^tV II . 2 ^ II 

L 1 ((0,T)xT) 



< ^ll^llL2( ( o,T)xr) ( / Afde ) " II < g £ > 117 



< C^fe -> 

as e — )• 0. Hence, according to (|4.6|) . we deduce that {£' e } £> o is uniformly bounded in 
(W{ ,w ((0,T)xr))'. □ 

As a direct consequence of Lemma 14. 1^ we have 

Lemma 4.2. E £ ^ E weakly in (Wq ,p )' , for some function E £ (Wq ' p )' mt/i p > 4, 
and (E, B) satisfies 

d t B + V x E = (4.7) 
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Proof. Indeed, the uniform bound on E £ in (Wq '°°) / and the Sobolev embedding 

W*' p ({0,T) xT)4 W^°°{(0,T) x T) 

for any 4 < p < oo imply that £7 e is uniformly bounded in (Wg' p ((0, T) x 7"))' and hence 
is weakly convergent in (Wq' p ((0,T) x T))' since (W 2 ' p ((0, T) x T))' with 4 < p < oo is a 
reflexive space. 
Next, since 

dB 
at 

holds in P'(R + x T), we take an arbitrarily test function </> G C^°(M + x M 3 ) to obtain 



/0 JT UL JO JT 

Hence, from (14.81). we obtain 



B e ■ -j^dxds + I I E £ -V x 4>dxds = 0. (4.£ 



* r [* f f f f d<h 

E ■ V x (ftdxds = lim / E £ -V x <j)dxds = lim / B e ■ —dxds 
JT ^°Jo JT £ ^°Jo JT dt 



dB 

B • —^dxds = — I I 6 • ——dxds 



(4.9) 



— *—dxds = — f 
10 JT dt ' Jo Jt^ dt 

Hence, from (|4.9p . we deduce that the limits (E,B) satisfy (|4.7p . □ 

Observe that, since £7 £ is convergent at least in the sense of distributions, we can 
conclude that £1 = 0. 



5. Vanishing of Conservation Defects 



Before stating the main result of the present section, we introduce a new class of bump 
functions as in [16]. For each C > 0, set 



T C = { 7 GT: IItV, < C} . 

Consider the transformation T defined by T7 = 1 — (1 — 7) 2 ; clearly T maps into T 2 c- 
Define 

f = TT 8 c T 16 , 
and notice that T / since Ts 7^ 0. For each 7 6 T, define 



j( z ) = 7 (z) + ( z _i)_Z. 



7 



Notice that 

supp7 C 

On the other hand, let 7 6 Tg be such that 7 = T7. One has 

c?7 



I 5 

2' 2 



3 5 
4' 4 



{I}- 



l-7(*) = (1-7) 



(l-7)-2(^-l) 



z > 



so that 



< 9(1-7), z>0. 



(5.1) 
(5.2) 
(5.3) 



(5.4) 
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Theorem 5.1 (Vanishing of conservation defects). Let'j £ T, and denote byr] = t/(£) 
any collision invariant (i.e. 7/(£) = 1 or r/(^) = £].,.•., £3 or else = \£\ 2 ) or any 

linear combination thereof. Then 



(5.5) 



X T) as e — > 0. 

Proof. We begin with the renormalized form (|2. 19[) of the Vlasov-Maxwell-Boltzmann 
equations $T2$) with T(z) = (z - 1)7(2) 



dt + h- Vx) (Mg b £ ) + e(eE E + £ x B £ ) ■ V C (M 5 ; 



c^7 



+ ee£ £ • £M<? e b - e (j(G e ) + (G e - l)^(G e )J # e • £F e 

= 1 / / (F' e f « - F £ F„) f 7 (G £ ) + (G £ - l)p-(G e )) bduM^. 
^ Js 2 Js. 3 V dz J 

Here, we used the decomposition f|3. 13|) . From (|5.6p . we deduce that 

d t + • (^) + ej {sE e + £ x B £ ) ■ V^Mg b )r,d£ 



(5.6) 



+ e / eE £ ■ ZMg b £V d£ - e / %E £ ■ £F £ r]d£ 

1 

= - <C q £ %rj >, 

where 

7e = 7(G £ ), 

and the function 7 is defined in terms of 7 by (j5.2[) 
Observing that 

(X x y) • Z = F • (Z x X) = X ■ (Y x Z), 

we have 



(5.7) 



/ 



(eE £ + ixB £ )-V i {Mg b £ ) 11 di 

eE £ ■ (v iV g b ) + J (£ x B £ ) ■ V iV g b Md£ 
eE £ ■ (Vt.rig b ) -B £ -Ux V^g b £ 



(5.8 



Notice that following the same line of the argument of Proposition 4.1 in |16j . it can be 
shown that 

1 







(5.9) 



^LOR+x-n 



as e — > 0. 

In order to estimate the L 1 -norm of the conservation defects, for the last two terms on 
the left-hand side of (|5.7p . we claim 



eE £ • (V 



0: 



(5.10) 
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eE £ ■ ^Mg b eV dC 



0; 



and 



%E £ ■ £F e ridZ - / E £ • CvMdC 
as e — > 0. Indeed, using the elementary bounds 

, , 27 

|7e|<9, |1 - Tel < 9, 0<G £ | 7e |<y, 
for the inequality (|5.10p . we have, 



< 



(jftfMdti) <C[ I {glfMdti 



(5.11) 



(5.12) 



(5.13) 



since 



f (V e r]) 2 MdZ < C 
Jr 3 



for all Tj 6 N(L) and where C is a positive constant. Hence, by the Cauchy-Schwarz 
inequality and the first statement in Lemma 13.21 one has 



eE £ -(V m 



< C 



e\E £ \ / {glfMdti 



-'loc\ 



\\g\ TP II 2 ll/i^ll 2 



< Cff2 -> 0, 

as e — )• 0. Similarly, for the inequality (|5.1ip . we have, 



< ( / m*Md£ >) { I {glfMdti ) < C ( / (<7°) 2 Md£ ) , 



since 



/ 

Jr- 



{trifMdt, < C 



for all r\ £ N{L) , where C is a positive constant. Hence, by the Cauchy-Schwarz inequality 
and the first statement in Lemma 13.21 one has 



eE £ ■ to&Mdt 



< C 



e\E £ \ U{glfMdA 



ill I 
<r rv? llc-9 p II 2 11,^11 2 

- ° fc H fc L?> (L 2 (T))" g e" L™(L 2 (Mdxd£)) 

l 



< C£2 -> 0, 



as e — )• 0. 
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It remains to deal with f|5. 12|) . To this end, we rewrite 



( %E £ ■ l;F £ nd( - I E e ■ £r/Md£ = ( {%- l)E £ ■ iF^di + e f E £ ■ £ We Md£ 

(5.14) 



h + h- 



Notice that from (|5.4p . we have 

|7 £ -1|<9(1-7(G E ))<9(1-7(G £ ))I 

for some 7 £ Tg and hence we can control I\ as, using F £ = MG £ , Lemma 13.21 and the 
fact < 1 - y(G B ) < 1, 



\h\\ L i oc( n + xT) < 9 



< 9 
+ 9e 

+ 9e 



E £ \\Cv\\l-l(G £ )\2G £ MdC 
E E \\^\\l-j(G e )\^Md^ 
|^ £ ||^||1-7(G E )|^ £ M^ 

Jv? 

\E e \\fri\g e Mdt 



^ oc ( K +x^) 



-xT) 



(5.15) 



^L( R +x^) 



(5.16) 



< i8^\\^E £ \\ LlciR+xT) (M 2 y \M\\\l UR+xT . L1{Md0) 

+ 9e\\E £ \\ L 2 oc{R+xT) (\£i]\ 2 ) 2 \\9e\\ L l c {R + xT^(Mdi)) 

< Cji + Ce -> 
as e — > 0. For 1%, we have 

II J 2||l1 oc (i; + xT) ^ V^\\V^E £ \\ L 2^ R+xT:i (|C?7|) 2 \\ge\\q oc {M.+ ;L 2 {Mdidx)) 

< Ce ->• 

as £ -)• 0. Adding ([BTTi]) . (IBTTBl and (IBTTBl) together gives (I5TT21) . Combining (l5?f|l - (l5TT2l) . 
the proof of (15.5f) is finished. □ 

Remark 5.1. According to Theorem 15. 1\ if 77 = 1 or 77 = | ^ | 2 , then the last term on the 
left hand side of (I5.5P will vanish; that is, 

E £ ■ (£> = E £ ■ (m 2 ) = 0, 

because 

(o = (m 2 ) = 0. 

This implies that the term E £ - < £77 > will only possibly appear in the conservation law 
of momentum. Hence, 

dt (gkk) + ■ (££fc<? £ 6 ) + eB £ ■ x V^ k g b e) ~ ae(E £ ) k -»• 0, (5.17) 
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in Lj oc (R+ x T) for all 1 < k < 3, since <^> = a = § (|£| 2 ). 

6. Proof of the Main Result: Theorem 12.11 

In this section, we will finish the proof of Theorem 12.11 via three steps. 

6.1. The Incompressibility and Boussinesq Relations. Let us start with considering 
the renormalized form of the first equation in ()2.2[) : 



ed t h £ + £ • V x h £ + ee (eE £ + £ x B £ ) ■ V^h £ - eeE £ ■ & = \^rQ{G £ , G e ), (6.1) 

N £ e z N £ 

where 

h e = -ln(l + ~eg e ] = -hiiV £ . 



e \ 3 / s 

Since h £ formally behaves like g £ for small e, it should be thought of as the normalized 
form of the fluctuations g £ . This means that, for every X € C l {T; L°°(Md£)) and every 
< s < t < oo, one has, 

e / < h £ (t)x > dx — e < h £ (s)x > dx — < h £ £ ■ V x x > dxdr 

Jt Jt Js Jt 

+ e [ [ e 2 E £ - < £h £ x > dxdr - e [ [ [ e(eE £ + £ x B e ) ■ V^xh £ Md^dxdr 

Js Jt Js Jt Jr 3 



Due to the fact 

and the entropy control 
one obtains 



(6.2) 

G K 



— < 3 

N £ ~ 

\\ £lE e\\Lf{L 2 {dx)) < C) 

J sE e • /e^\ dxdr -+ 0, 

as e — t- 0. On the other hand, since as stated in the last statement of Theorem 13.11 
(cf. also Corollary 3.2 in [2]) that h £ has the same limit g as the sequence g £ in w — 
Lf oc (dt; w 2 L (M d^dx)) , one deduces that 



and 



e J ( h e(t)X) dx-e (h £ (s)x) dx -> 0; 
J J^e 2 E £ -(Zh £X )dxdT^0; 

[[[ e{eE £ + £ x £ £ ) ■ V '^xKM d^dxdr -»• 0, 

Jt J*. 3 
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as £ — > 0, thanks to the uniform bounds 

i 

II £ 2 ^11 L- + (L2(T)) < C > \\ B e\\L^ (L 2 (T)) < C. 

Taking the limit in (|6.2[) as e tends to zero while using Theorem 13. II to establish the limits 
of the terms involving h £ and q e respectively yields 

{g£ ■ VxX> dxdr = ((qx)) dxdr; 

t Js Jr 

hence, the limiting form of (16. ip is 

i-V x g = J J qb(Z t -Z,u)dbjM t (% t . (6.3) 

Since q is in L 2 (dAidx), then for every r\ = r/(£) in L 2 (dA4), an application of the 
Cauchy-Schwarz inequality shows that ((rjq)) is in L 2 {dx). By a repeated application of 
the <i.A/f-symmetries in Theorem 13. H one has that, for any r\ in L 2 (dM), 

({vq)) = lUv + v*-v'*-ri > )q}- (6-4) 

Successively apply the identity (|6.4|) for 77 = 1, £, ||£| 2 and use the microscopic conservation 
laws (|1.2p to obtain 

((g)) = o, «fr» = o, (^l 2 ff) = o. 

Since these 77 are also in L 2 (Md^), it then follows from the limiting Vlasov-Maxwell- 
Boltzmann equation (|6.3p that g satisfies the local conservation laws of mass, momentum, 
and energy: 

div x (£g) = 0, div x (£ ® £g) = 0, div, (d\£\ 2 g) = 0. (6.5) 



, '2 

Theorem 13.11 states that 5 has the form of the infinitesimal Maxwellian 

g = h + u-t + 9(~\£\ 2 -~ 

Substituting this into (16. 5p . the local mass and energy conservation laws yield the incom- 
pressibility relation for the velocity field u while that of momentum yields the Boussinesq 
relation between h and 9: 

div x u = 0, V x (h + 9) = 0. 

6.2. Proof of Convergence to Incompressible Electron-Magnetohydrodynamic- 
Fourier Equations. Throughout this subsection, it is assumed that the bump function 
7 belongs to T (defined by (|5,ip ). Using Theorem 15.14 the classical Sobolev embedding 
theorems, and the continuity of pseudo-differential operators of order on W s ' p for 1 < 
p < 00, one sees that, for all s > 0, 

a t p + pv* ■ \ ( U ® e - ^iei 2 /) gl 

(6.6) 



+ eP[B £ -{CxV^g° e )) -aePE e 
► 
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* L L(dt;W-J A (M. 3 )), and 



s,l 

l.,.o A A _ 1 /.A 



ft ^ ^lel 2 - 1 J flg^ + V. ■ - ^ (^| 2 - lj fl^ (6.7) 

in Lj oc (dtdx) as e — > 0. Here, the operator F is the Leray projection, i.e. the L 2 (dx)- 
orthogonal projection on the space of divergence-free vector fields. In (|6.7p . we used 

^xv c Qiei 2 -i) =o. 

By Theorem 13.11 and Proposition 13.21 pick any sequence e n — > such that 

g b £n ^g in w* - L? {L 2 {M d&x)) , (6.8) 

7 Bn «L-^9 in ^-^(^(^^((l + lel 2 )^))- (6.9) 

In this section, we deal exclusively with such extracted sequences, drop the index n and 
abuse the notations g £ , g £ , g £ , q £ and so on to designate the subsequences g £n , g £ , g £ , 
q £n . Set u and 6 the limiting fluctuations of velocity and temperature fields defined by 



tg b £ )^u, in w*-L?(L$); (6.10) 

l\Z\ 2 - Agty^O, m w*-L?{Ll). (6.11) 

The second entropy control in Proposition 13.21 implies that g £ and g £ have the same 
limit g in w — Lj oc (dtdx; L 1 (M<i^)); hence the Boussinesq relation and the incompressibility 
condition hold: 

div^u = 0, 9+<g>= 0. (6.12) 

Denote by ? either the tensor <3> or the vector Since L is self-adjoint on L 2 (Md£) so 
that 

\ ((L,)gl) = \ (,{Lgt)) = \ ((^ + £ - £ - £'))) 



\(g b e + al - g h l - g h J) + (g b £ gl - g h Jg h J) 

b 



(6.13) 



+ (sQ{ge,g 



The first term on the last right hand side of (|6.13p converges to the diffusion term while the 
second term converges to the convection term in the incompressible system (jl.4p . These 
limits are analyzed in the next two lemmas. The convergence to the diffusion term is 
obtained by an argument that closely follows |16| . except that the present work should 
pay additional attention to the Maxwell effect. This apparently minor difference makes 
our analysis slightly more difficult than that in |16j . 

Lemma 6.1. Define 



Then, as e — > 0, 



— ($:£$>, « = — (6.14) 



\ ((mg b e) - (*Q(gl9t)) -mv,u + (V,u) 
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in w — Lj oc (dtdx). 

The convection term is the nonlinear part of the limiting system and its convergence is 
therefore the most difficult to establish. The analysis below rests not only on all a priori 
estimates and the arguments in [16], but also the compactness of the moment of g\ in £ 
which is stated in Lemma 16.31 below. 

Lemma 6.2. The following convergence hold in the sense of distributions on M + x M 3 ; 

PV* ■ ($Q(g b ,g b )) -> PV X ■ (u ® u), 

as e — > 0. 



6.3. The Lorentz Force Term. The key result of this subsection is to deal with the 
convergence of the Lorentz force term. To this end, we first state the following compactness 
about the moment of g £ in £. 

Lemma 6.3. Let 7 € T be the same as in (|3.12p and the hypothesis (H2) hold. Then, the 
family g\ has the following property: for each sequence e n — > 0, each function \ = x(0 
as |£| —> 00, each T > 0, there exists a function rj : 1— > R+ such 
- and 

J \{9 b s n X) (t,x + y) - (g b £nX ) (t,x)\ 2 dxdt < V (\y\) 
for each t/6l 3 such that \y\ < 1, uniformly in n. 

Proof. For any 7 G Y, since F e is a renormalized solution of (|2.4p relatively to M, using 
the nonlinear function T(z) = (z-l)j(z) in the renormalized formulation (|2.19p . we obtain 



such that ■ 
that lim 2 ^ + Vi z 



(ed t + Z-V x )g b = [ [ q e %bduM*d£, - ediv f (e(eE £ + f x B £ )g b \ 
Jm.3 Js 2 v ' 

+ e%eE £ -£G £ , 

with 7 defined in terms of the truncation 7 by (|5.2p . Denoting 

7, if |/|<£; 

/ A L = I L, if / > L; 

-L, if f<-L 

for every L > 1, we deduce from (16.15P that 

(e$ + £ • V s )(s* A L) = Qf ^ ^ q^bduM^^j l { \ g b\< L] 

- ediv f (e(e£ 6 +(xB £ ) (<? e 6 A 
+ e7 £ e J B £ -£G e l { | 9 _6|< L} , 



(6.15) 



(6.16) 
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Furthermore, for every N > 1, we decompose g h e A L as 



g b e AL = gb + gb, gf = 0, 



with 



and 



(ed t + £ ■ V x )g\ = N J q £ %bduM^?j l{| ff ||<L}l{|^[>iV}, (6-17) 



where 



(ed t + £ • V 

- ediv € (e(eE £ + £ x B £ ) [g b £ A 
+ e%eE e ■ £G e l { \ g b\< L} , 



A e = I / q £ %bdwM*d£*. 
Is 2 



(6.18) 



Step 1: Control of g\. From (16.17p . if we denote 

S e = (J J q £ %bdu)M*d£^J l{\ Ae \>N}l{\ g b\<L}, 

then we obtain 

gl(t,x,0= [ e S e (t-es,x-st,Ods. (6.19) 
J o 

Notice that, since |7 e | < 9 and q e is weakly compact in ^(dtdxdAi), S £ is uniformly 
bounded in L l (dtdxMd£). Therefore, 

Step 2: Compactness of g\. Setting 

Se = {LL qe%bdwM * d v i{iA ^ N}i {\ 9 i\<L} 

- ediv f [e(eE £ + £ x B £ ) [g b e A l)) 

Notice that |7 £ G £ | < and hence, by the interpolation between L 1 and L°°, we have 
(/a /a 1^ bdujM * d ^ 1 {l^l< Ar } 1 || 9 b|<L} + e ^ eEe ' ^ e L 2 (dtdxMdx) 

and 

div e (e(££ e + (x B 6 ) Ail) 6 L 2 (dtdx; H' 1 (d£)) . 
Thus, from (|6.18p . we obtain 

(edt + £ • V x )^ = 4 € L 2 (dtdxMdx) + L 2 (dtdx; H- 1 (d(). (6.21) 



26 



XIANPENG HU AND DEHUA WANG 



Applying the averaging theorem in JSj [2], we deduce from (|6.2ip that, for all x(0 such 



that 



x(0 
T+HF 



as |f | 



oo, 



9sX 



i < Cn l, 
L 2 (0,T;H?(T)) 



(6.22) 



where Cn,l depends only on N,L. This yields the compactness of < g b e x > in space; 
namely, there exists a function n : R + i— > R + such that lim 2 _ s> o+ v{ z ) = 



L 2 ((o,T)xT) 



<^(M)- 



(6.23) 



Step 3: Compactness of <^ A L. From (|6.2ip and the weak compactness of q £ in 
L 1 (dtdxdAi) , we have, for large enough N, g b (t,x,^) 

like. Thus, by (|6.23p . there exists a function i] : M + i— >• K 
^AL) X )(t,- + y)-([g b £ AL) X ) (i, •) 



can be as small as we 

^(dtdxMdt) 
l + with lim z _ i> o+ 77(2) 



0, 



L!((0,T)xT) 



<^(M)- 



(6.24) 



Using the hypothesis that \ (de) 2 \ is relatively compact in w — L 1 (dt(l + \£\ 2 )Md£dx), 



! "I 

we deduce easily that there exists a function rj : i— )■ K+ with lim z _> + r /( z ) = 0> 



^Al) X )(i,-+y)-((^AL) X )(V) 



L 2 ((0,T)xT) 



<^(M)- 



(6.25) 



Step 4: Compactness of g^. Due to the hypothesis that {(^e) 2 }^ y * s relatively 

compact inw - L l {dt{\ + \^\ 2 )Md^dx), for every (3 > 0, there exists an integer L > 1 such 
that 

<C/3, 



](s£al) X ) (*,-)-( (flg) x) (V) 
uniformly in e. Thus, for such /3 and L, we have 



L 2 ((0,T)xT) 



s^Ai )x){t,- + y)-(g b £ x)(t,- + y) 



and 



L 2 ((0,T)xT) 

< cp 



<cp, 



Ai)x)(V)-(^)(V)| L2((oT)xr) 

uniformly in e. Hence, the above two inequalities, combining together with (|6,25p . imply 
there exists a function n : R + 1— > M + such that lim z _>o+ r/(z) = and 
rT t 2 

(g h £nX ) (t,x + y) - (g b £nX ) (t,x)\ dxdt < V (\y\) 



JT 



for each !/£l 3 such that \y\ < 1, uniformly in n. 

Now, we are ready to prove the convergence of the term of Lorentz force. 
Lemma 6.4. The following convergence holds in the sense of distributions on 



□ 



P 3. 



5 £ -(?x V^ k g° £ ) -> (B x (V x fl)) fc , 

as e — > 0, for all 1 < k < 3. T/ie notation a k stands for the i-th component of the vector 
a. Further, we have j = eu. 
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Proof. For any 1 < k < 3, = e^, where {ek}k=i ^ s the standard basis for M 3 . This 

implies 

x V^) = (^) x efc . 
Then, we can rewrite B £ - < £ x V^k9 b £ > as 



Defining 



<e(l+eflg)> 



(6.26) 



Je 



since (0 = 0. Then, we have 



Je e 



-»■ 0, (6.27) 

Lf (HftfcATde)) 

as e — > 0. Indeed, from the definition of g £ , we know that eg c e is uniformly bounded in 
L ( ^°(L 2 (dxMd£ i )) while from the second statement of Lemma [3.21 g £ is uniformly bounded 
in L^(L 1 (dxMd£ t )). Thus, by the interpolation between 1? and L , we deduce that 

11 ^ £U L?°(L?(dxMd£)) — ' 

for some constant C > 0. Therefore, we have 

~ 9ei 



n b - il 

Je e 



Lf > {L 1 {dxMd£,)) 



Lf (L 1 (dxAfd$)) 



£f (^(tteAfdC)) 



< \\Z9 



£ "L^°{L 1 2(dxMd£,)) 



<IC| 3 > = 



< C£2 



0. 



as e — > 0. Hence, (|6.27p . combining with the weak convergence of m L^°(L 2 (dxMd^)) 

and the uniform bound of {j b } in L^°(L 2 (dxMd^)), implies that j h £ converges weakly to j 
in Lf(L 2 (dxMd£)). Note that if = e(g E £), we have j = eu. 
Notice that, (|2.2c[) implies 

d t B £ = -V X E e e L°°(0,T;W-^ 2 (T)) C L 1 (0, T; W~ S,1 (T)), 

for some s > 4 large enough, and is bounded in -L 1 (0, T; W~ S ' 1 (T)) uniformly in e. On 
the other hand, Lemma 16.31 with x(0 = 6c implies that for each T > 0, 

J J I (*, » + y) x e fc - (^) (t, x) x e fc | 2 dx(it < rj(\y\), 



for each y G IR 3 such that |y| < 1, uniformly in e, where rj is a function . 
lim z ^o+ r /( z ) = 0- Hence, by Lemma 5.1 in [26], one has 



(6.28) 
satisfying 



(B £ ) • x e fc J -> 5 • (eu x e fc ) = B ■ (j x e k ) = (B x j)fc, 

in the sense of distributions. The proof is complete. 



(6.29) 
□ 
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